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(Abstract) Recently Pao et al proposed a method to solve the transient response of truss type structures under impact load: the 
method of reverberation ray matrix (MRRM). This is a novel matrix method in frequency domain that takes into account the 
multi-scattering of axial and bending waves at the joints and derives scattering coefficients upon equilibrium and compatibility 
conditions of the joints. The present work develops formulation of MRRM in terms of Euler-Bernoulli beam theory and analyzes the 
transient response of the structure subjected to step or impulse displacement load. The efficiency and accuracy of MRRM is verified 
by comparison of results of MRRM and Mindlin' s solution for a single beam, and that of FEM for the transient response of axial and 
shear forces, bending moment, displacement, velocity and acceleration of a 17-bar rigidly jointed truss. In general, it is found that 
the result of MRRM still coincides with that of FEM. While, the solving process of MRRM is more simple and quick than that by 
FEM, and the MRRM needs less computing elements. 



1. INTRODUCTION 

In recent years, there has been an increasing interest on the 
dynamic response of truss or frame structures as need arises in 
outer space engineering, for instance, for the purpose of 
controlling local disturbance initiated by incident impact or 
monitoring of structure integrity of lattice type space station, as 
well as in civil engineering. There are mainly two approaches 
for the analysis of dynamic response of truss or frame. One is 
vibration analysis based on mode superposition approach, 
which finds the vibration modes of the structure and express the 
dynamic response by superposition of all components in the 
modal coordinates. The other one is to study the dynamic 
response of a structure by elastic wave theory. If the response is 
induced by impact load, then the initial stage of the transient 
response contains many high frequency components. The 
analyses based on wave propagation approach are shown to be 
superior to the other methods. 

One of the early research work with ray tracing method was 
reported by Boley and Chao' 1 '. They studied how stress waves 
are generated by impact force at a joint of a truss, propagate 
along the members emanating from the loaded joint, and scatter 
from the joints at opposite ends. It was noted that by 
considering all possible paths of propagation the dynamic 



response of a truss can be described as a superposition of waves 
that reverberate within the truss. Only the first two or three 
reverberations of axial waves in their example were calculated 
due to lack of computer at that time. The scattering of incident 
axial wave at a rigid joint of three bars were analyzed and 
experimentally investigated by Desmond . Experimental 
methodology for analyzing the effect of structural joints on 
wave propagation was presented by Doyle and Kamle 131 . A 
matrix formulation for planar frame was proposed by Nagem 
and Williams 141 which combines the transfer matrices with joint 
coupling matrices to form the system matrix that characterizes 
the steady state response of the structure. With the system 
matrix they have calculated the natural frequencies of a sample 
structure. 

Recently, Howard and Pao' 5 'proposed a novel matrix method, 
i.e. the MRRM, which investigates the dynamic response of the 
planar trusses in terms of axial waves propagating along 
structural members and scattering at the joints. In this method 
the scattering coefficients representing the reflection and 
transmission of axial waves at each joint are derived from the 
equilibrium and compatibility conditions at the joint. 
Comparison between theoretical results and experiments 
indicate that the axial wave theory is only valid for the response 
at the very early time. In a following paper, based on 
Timoshenko beam theory the method is modified by Pao, Keh 
and Howard' 6 ' to include the effect of flexural waves and mode 



conversions between axial waves and flexural waves at the joint. 
Comparison to experimental data of the model truss under a 
step load shows good agreement for the early as well as 
considerably long time responses. Afterwards, Pao and Sun' 7 ' 
applies MRRM to analyze the dynamic behavior of pin-jointed 
trusses subjected to suddenly applied force. It is found that the 
amplitude of dynamic bending strains in the members of 
pin-jointed trusses, as in the rigid jointed trusses, are 
comparable with the dynamic axial strains in the same member. 
Such a large dynamic bending strain in a truss member, which 
is contradictory to the static behavior of the truss, could only be 
caused by the inertial force of the mass in structural members of 
the truss. 

In Ref. Pao et aF ] compared the method of transfer matrix 
(MTM) with MRRM in detail, concluded that MRRM is a 
viable alternative to the solutions of initial value and two point 
boundary problems of unidirectional wave motions and could 
be extended further to analyze wave propagation in two 
dimensional structures. Additionally, the MRRM is employed 
to obtain frequency response function (FRF) of displacement of 
a frame under the action of a unit impulse load by Miao et al '. 
The natural frequencies of the frame are determined from the 
FRF, the vibration mode is retrieved from the adjoint matrix of 
the coefficient matrix of the governing equations of MRRM. 
The MRRM has advantage over numerical methods, such as 
FEM. 

The investigations of dynamic response of truss type 
structures with MRRM are still limited in the literature. 
Nothing else than the responses of bending and axial strains of 
trusses subjected to Heaviside step force have been reported. 
However, some of the engineering problems, such as sudden 
subsidence of foundation of engineering structures, can be 
categorized as suddenly applied displacement at support or 
joint of the structure. In the present study, the formulation of 
MRRM is re-derived based on Euler-Bernoulli beam 
theory that is a reduced case of the truss made of 
laminated beams. The latter is to be analyzed in our 
further research. The transient response of planar rigidly 
jointed truss subjected to step or impulse type 
displacement load is addressed with MRRM. The 
efficiency and accuracy of MRRM are demonstrated by 
comparing solutions of MRRM with results obtained by 
FEM. 

2. THEORETICAL ANALYSIS 

A planar truss-type structure is an assemblage of rods or 
beams pin-connected or rigidly connected at joints. In the 
following only rigidly jointed truss is considered. We define 
joints by letters (I ,J ,K ,...) or numbers ( 1, 2, 3, . . . )and 
members by two letters or two numbers ( IJ , JK 12, 23, 
. . . ). As shown in Fig.l, member JK means one with joint J 
as starting point and joint K as the end. The origin of a local 

JK JK JK 

coordinates x — y locates at J with x -axis laid along 
JK. The number of members meeting at the joint J is n J .X-Y 
is the global coordinate system. 



The governing differential equations based on 
Euler-Bernoulli beam theory for a uniform beam in frequency 
domain can be written as: 



E d 2 u(x,co) 
u(x, co) H — — = 



pco 



dx 1 



EI d v(x, co) 

v(x,G)) —j —j = 

pAco dx 



(1) 
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where, E , p , A , / are elastic modulus, density, cross 
sectional area and moment of inertia; u , v are axial and 
transverse displacements in the local x-y coordinate system. 
Then, axial force F , shear force Q bending moment M and 
angle of rotation cp of the member can be expressed as: 

du d\ d 2 v dv 

F = EA — , Q = -EI — -, M = EI — -, <p = — ; 

dx dx dx dx 



The solutions for u(x, CO) and v(x, CO) are assumed as: 



u = a x e 1 + d x e 



iK^x , i —ifc-yX , ik-iX , i —ik-iX 

v = a 2 e 2 +a 2 e 2 +a 3 e +a 3 e 



(4) 
(5) 



dv 



<p = — = ik 2 a 2 e' kl ' - m. 2 u 2 
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where a, (co) , d i (co) (i =1,2,3) are the amplitude of 
arriving wave and departing wave and kj(co) (i =1,2,3) the 
wave number of the i-th mode, respectively. Substituting 
eqs.(4)~(6) into eqs.(l)~ (2), the wave number can be written 
as: 

(7) 



The external forces and moment applied at the joint J are 
denoted as f x , f Y , rri . In the global coordinate system, the 
displacement in the X and the Y directions and the rotation of 
the joint J are denoted as U J (co) , V J (co) and <£> J (co) , 
respectively. 6 is the orientation of the member JK with 
respect to X axis. The equations for the balance of forces along 
the X and the Y directions and the moments at the joint J are: 

]T [F M (0, co) cos 6 JK - Q JK (0, co) sin 6 JK }+ft=Q (») 



]T [F JK (0, co) sin 9 JK + Q JK (0, co) cos JK ] +// =0 (9) 



^M J *(0,») + m J =0 (io) 

K=l 
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The compatibility conditions of the displacements along the X 
and the Y directions and the rotation at the joint J are: 

u JK (0,co) cos JK -v JK (Q,co)smO JK =U J (co) 

(x=i,2,...y ) (ii) 

u JK (0, co) sin 9 JK + v M (0, co) cos JK = V J (co) 

(K=\X.-.,n s ) (12) 

cp M (0, CO) = J (CO) (K = l2,...,n) (13) 

The equilibrium equations and displacement compatibility 
conditions at joint J of the truss can be expressed in the matrix 
form as Ref.[7]: 

,i J fd j i rA,i\ ,i r oi 

(14) 
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d = [I-R] *s 



(18) 



where R = SPU is named as the reverberation ray matrix for the 
structure. It is noticed that the Fourier inverse transform of the 
matrix [I — R] 1 involves an infinite number of poles along the 
real axis of the complex CO plane. To circumvent this difficulty, 
Pao et al [b] proposed that [I — R] 1 is expanded into a 
Neumann series. 

[I-Rr 1 =I + R + R 2 +--- + R"+--- 

(ly) 

Substituting Esq. (18) and (19) into Esq. (3) ~ (6), we can 
obtain the forces, the moment, the displacements and the 
rotation at any location of the member JK in the frequency 
domain. The unknowns of the joint J {H J } can also be solved, 
since d and a have already known. Then, applying Fourier 
inverse transform, we can obtain the dynamic response of 
above variables in the time domain. 



For the above equations, vectors and matrixes are in bold font. 
They are all functions of frequency co. {dJ} and { aJ} are vectors 
of amplitudes of departing and arriving waves for all members 
jointed at the joint J. The superscript J appearing in matrix [D] 
and [A] means they are the coefficient matrix of the joint J. 
Those coefficients inside square brackets consist of material 
parameters and geometric parameters of the truss. {HJ} and {fj} 
consist of the unknowns and known quantities of the joint J in 
the global coordinate system respectively. For non-constrained 
degree of freedom of the joint, the displacement (angle of 
rotation) is unknown while the force (moment) is given. For 
constrained degree of freedom of the joint, the constraint force 
is unknown while the displacement of the joint is given. 
Solving eq.(14), the amplitudes of departing waves can be 
expressed by th^ amplitudes of arriving waves and the local 
source vector s at the joint J . It can be written in the matrix 
form as: 

d J =SV+S J (i 5 ) 

where S J is named as the local scattering matrix of the joint J . 
Assembling amplitudes of departing waves for all joints into a 
global vector d , and amplitudes of the arriving waves into a , 
the global equation relating a and d can be written as 

d = Sa + S (16) 

The matrix S is referred to as the global scattering matrix, and 
the vector s as the global source vector. It is shown in Ref. [6] 
that the vector a is related to the vector d as 

a = PUd (17) 

in which P is the propagation matrix and U is the permutation 
matrix. Equation (16) then becomes 

d = SPUd + s 

This equation can be reduced to 



3. EXAMPLE OF A SIMPLY SUPPORT BEAM 

To verify the theory and the computer code of MRRM, we 
first analyze a simple structure by the established theory. A 
simply supported beam is shown in Fig2. Assume the beam is 
made of aluminum alloy with square cross section. The side 
length of the cross section is 0.635cm. The elastic modulus of 
aluminum is 69GPa, density is 2700Kg/m 3 , and axial wave 
speed is c 1 =5039m/sec. Take l o =0. 1016m as a unit of length, 
the length of the beam is / =0.4064m (/=4/ ). A step 
displacement expressed using Heaviside function S-H(t) 
((5=0.0 l/o) is applied to the right end of the beam. The unit of 
time is taken as t = l /c l = 20.16//S , which is the time 
needed for axial wave to propagate one unit of length. With 
MRRM, the beam is discretized into two equal-size elements. 
Joint 2 is at the mid-span of the beam. 

Applying Mindlin's method' 101 , the transverse 
displacement of beam has the form 171 : 



x,0 = [7 + V (-1)" — sin^cos((y„f)]^ (20) 
/ t~7 nn I 



where CO n = (tlTT I 1) 2 ^EI I pA . The results of transverse 
displacement at mid-span of the beam obtained with MRRM 
are compared with that of eq.(20) in Fig. 3. In the figure the 
solid curve v23 is the result of displacement at the starting joint 
2 of member 23 in the local coordinates, obtained from eq.(5) 
and evaluated at x = 0. The dashed curve V2 is the result of 
vertical displacement of Joint 2 in the global coordinate system, 
which is obtained as joint unknown in { H J } . V2 should be equal 
to v23(0) by the compatibility condition eq.(12). The dotted 
curve is the result of eq.(20) with x =0.5/ and 40 terms are 
included in series expansion. It is found that the first two curves 
resulted from the two solutions of MRRM are superposed and 
coincide very well with that of Mindlin's method. 



4. TRANSIENT RESPONSE OF A 17-BAR 
FRAME SUBJECTED TO SEMI-SINE 
DISPLACEMENT IMPLUSE 

In this section the laboratory model of the 17-bar aluminum 
rigidly-jointed truss, which is experimentally and analytically 
investigated in Ref.[7], is used as example. The material 
properties and cross sectional dimension are the same as the 
beam considered in the previous section. As shown in Fig. 4, the 
length of the horizontal member (bar 35, 24, etc.), the vertical 
member (bar 34, 56 etc.) and the diagonal member (bar 45 etc.) 
are3/ , 4/ and5/ , respectively. The total span of this frame 
is L =4x3/ =1.219m. The unit of time is also t Q . As shown in 
Fig. 5, the impact semi-sine displacement applied at joint 6 is 
V 6 = 8 ■ sin(^//r) , in which 8 = 0.0 1/ , T = 9t . From 9 
to 30 units of time, displacement amplitude is zero. In the 
frequency domain the semi-sine displacement impulse can be 
expressed as 

W= ; S '* /T , (1 + ^) (2D 



(nix) 1 



00 



According to the above theory, we can calculate the transient 
displacement, force and moment at any location of each 
member, as well as the joint unknowns by inverse fast Fourier 
transform. 



4.1. EQUILIBRIUM OF FORCES AND BENDING 
MOMENTS AT A JOINT 

In the general theory the forces and the bending moments of 
all members connected at a joint are balanced at all times, as 
stated in eqs.(8)~(10). In Fig. 6, theXcomponent of the end 
force at the stating point of members 53, 54, 56 and 57 (i.e. at 
the point connecting to joint 5) are plotted. They are denoted as 
F x 53, F x 54, F x 56 and F x 57 respectively in Fig.6. 
Each of them is the sum of projections of axial force and shear 
force for the corresponding member when x = in the local 
coordinates, as shown by the two terms in bracket of eq.(8). The 
Y component of the end force of these members are denoted as 
F Y 53 , F Y 54 , F Y 56 and F Y 57 , and plotted in Fig.7. The 
unit of force is f Q = 712N . In Fig. 8 the bending moments of 
the four members at joint 5 are plotted. The unit of bending 
moment is m = 72.36N ■ m . The summations of these 
forces inXor in Y directions and the moments should be zero 
according to equilibrium condition of eqs.(8-10). It is found 
that the sums of the forces in either X or Y direction as well as 
the bending moments from the four members connected at joint 
5 are balanced very well, which are depicted by the dashed 
straight line along the time axis. 

4.2. COMPARISON WITH THE RESULTS OF FEM 

In this section the results of MRRM are compared with that 
of FEM. In FEM, each member is discretized into 100 elements. 
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While for MRRM the variables at any position of the 
member can be evaluated without discretization of any member. 
Figures 9-11 show the comparison of axial force, shear force 
and bending moment at the middle point of member 57 
obtained from MRRM and from FEM. There are 2 curves in 
each figure. The curve in solid line is the result of MRRM, that 
in dotted line is the result of FEM, the unit of force and moment 
are also /„ =712N and m =72.4N-m . The results of 
MRRM and FEM show good agreement. 

Figure 12 show the comparison of vertical displacement at 
joint 5 resulted from MRRM and FEM. In the figure v53, v54, 
v56 and v57 are the vertical component of displacement of the 
corresponding member connected at joint 5. Each of them is the 
sum of projections in Y direction of axial and transverse 
displacement u and v for the corresponding member when x = 
in the local coordinates. They are found to be dispersed slightly. 
V5 is the vertical displacement of joint 5 calculated as joint 
unknown variable from {H J }. K5_FEM is the curve of vertical 
displacement of joint 5 obtained from FEM. In general these 
results coincide very well. It is noticed that the maximum 
displacement of joint 5 along the Y direction is about 0.0 1 2/ , 
which is larger than the amplitude of displacement inputted at 
joint 6 (0.01 /„), 

Shown in Fig. 13 are vertical velocities calculated from 
MRRM and FEM. One unit of velocity is v = 5039m/s . The 
first four curves denoted as ve/53, ve/54, ve/56 and ve/57 are 
corresponding responses in Y direction of the global coordinate 
system of the four members at joint 5, which are converted 
from solutions of axial and transverse velocities in the local 
coordinate system. The curve denoted as ( / e/5_FEM is the 
response of vertical velocity of joint 5 resulted from FEM. It is 
found that there is no response during the initial 4t Q time 
before substantial response is observed. 4t Q is the time for the 
disturbance to travel from joint 6 to joint 5. 

Figure 14 compares results of vertical acceleration at joint 5 
obtained from MRRM and FEM. One unit of acceleration is 
a = 2.49xl0 8 m/s 2 . The first four curves denoted as 
acc53, acc54, acc56 and acc57 are corresponding accelerations 
in Y direction of the global coordinate system of the four 
members at joint 5, which are converted from solutions in local 
coordinate system. It is found that they are completely 
superposed. The curve denoted as ^4cc5_FEM is the result of 
FEM. In general, it still coincides with the results of MRRM. 
While, it also can be found that the solving process of MRRM 
is more simple and quick than that by FEM, and the MRRM 
needs less computing elements. 

5. TRANSIENT RESPONSE OF THE 17-BAR 
TRUSS SUBJECTED TO STEP 
DISPLACEMENT 

In figures 15 and 16 the vertical velocity and vertical 
acceleration at joint 5 are presented for the 17-bar rigidly 
jointed truss when it is subjected to a step vertical displacement 
8 ■ H(t) (<5=0.01/o) at joint 6. It is found that in these two 
figures the first four curves denoted as ve/53, ve/54, ve/56, 



ve/57 in Fig. 15 and acc53, acc54, acc56, acc51 in Fig. 16 
respectively, which come from the solution in the local 
coordinates of MRRM of each member connected at joint 5, are 
completely superposed. However, more disagreement between 
solutions of MRRM and FEM for the truss subjected to step 
displacement are observed than the truss subjected to semi-sine 
displacement load as presented in the previous section. 

6. CONCLUSION 

The MRRM, which was recently proposed by Pao et af \ has 
been proved both theoretically and experimentally an efficient 
approach for transient dynamic analysis of truss type structures. 
This is a matrix analysis method in frequency domain that is 
successful in resolving singularity problem by expanding the 
singular matrix in the left side of eq.(19) into Neumann series 
before variables are inversely transformed by IFFT to time 
domain. In the present work the MRRM is substantiated in 
dynamic analysis of rigidly jointed truss where the load is given 
by a step or an impulse of displacement at a joint. The analysis 
and the computing of MRRM are verified first by the 
agreement between the results of MRRM and Mindlin's 
method of the simply supported beam subjected to end step 
displacement. Then the balance of the forces and moments at 
the ends of members connecting at a joint of the 17-bar 
truss are demonstrated. The analysis and computing are 
further verified by agreement between the results of 
MRRM and that of FEM for force, moment, displacement, 
velocity and acceleration of the 17-bar truss subjected to 
step or semi-sine displacement load. Compared with the 
FEM, the same results can be obtained by MRRM with 
less computing elements, and the solving processes are 
more simple and quick. 
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Figure captions 

Figure 1 Local and global coordinate systems of the truss 
Figure 2 Simply supported beam subjected to step 
displacement at the right end 

Figure 3 Comparison of vertical displacement at the middle 
point of the beam by MRRM and by Mindlin's method 
Figure 4 17-bar aluminum truss subjected to displacement load 
at joint 6 

Figure 5 Load history of semi-sine displacement 
Figure 6 Equilibrium of forces in X direction at joint 5 
Figure 7 Equilibrium of forces in Y direction at joint 5 
Figure 8 Equilibrium of bending moments at joint 5 
Figure 9 The dynamic response of axial force at the middle 
point of member57 

Figure 10 The dynamic response of shear force at the middle 
point of member57 

Figure 1 1 The dynamic response of bending moment at the 

middle point of member57 

Figure 12 Vertical displacement at joint 5 

Figure 13 Vertical velocity at joint 5 

Figure 14 Vertical acceleration at joint 5 

Figure 15 Vertical velocity at joint 5 when the truss subjected 

to step displacement 

Figure 16 Vertical acceleration at joint 5 when the truss 
subjected to step displacement ^ 
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Fig.l Local and global coordinate systems of tlp-tijtjs^ 



c 



HO 



1/2 



I 



Fig. 2 Simply supported beam subjected to step displacement at 
the right end 
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Fig. 3 Comparison of vertical displacement at the middle point 
of the beam by MRRM and by Mindlin' s method 




Fig. 6 Equilibrium of forces in X direction at joint 5 
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Fig.4 17-bar aluminum truss subjected to displacement load at 
joint 6 
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Fig. 7 Equilibrium of forces in Y direction at joint 5 
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Fig. 5 Load history of semi-sine displacement 
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Fig. 8 Equilibrium of bending moments at joint 5 
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Fig. 9 The dynamic response of axial force at the middle point 
of member57 
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Fig. 10 The dynamic response of shear force at the middle point 
of member57 
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Fig. 1 1 The dynamic response of bending moment at the middle 
point of member57 





v53 




v54 




v56 




v57 




V5 




V5_FEM 











10 15 
^0 



Fig. 12 Vertical displacement at joint 5 
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Fig. 13 Vertical velocity at joint 5 
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Fig. 14 Vertical acceleration at joint 5 




Fig. 15 Vertical velocity at joint 5 when the truss subjected to 
step displacement 




Fig. 16 Vertical acceleration at joint 5 when the truss subjected 
to step displacement 



